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We implement a technique to characterize the electromagnetic properties at frequen¬ 
cies 100 to 165 GHz (3 cm“^ to 4.95 cm“^) of oriented smectite samples using an 
open cavity resonator connected to a sub-millimeter wave VNA (Vector Network An¬ 
alyzer). We measured dielectric constants perpendicular to the bedding plane on 
oriented Na"*" and Ca’''’''-ion stabilized smectite samples deposited on a glass slide at 
ambient laboratory conditions (room temperature and room light). The clay layer 
is much thinner (~ 30 /rm) than the glass substrate (~ 2.18 mm). The real part of 
dielectric constant, Cre, is essentially constant over this frequency range but is larger 
in Na"*"- than in Ca’''’''-ion infused clay. The total electrical conductivity (associated 
with the imaginary part of dielectric constant, e^m) of both samples increases mono- 
tonically at lower frequencies (< 110 GHz), but shows rapid increase for Na^ ions in 
the regime > 110 GHz. The dispersion of the samples display a dependence on the 
ionic strength in the clay interlayers, i.e., ^-potential in the Stern layers. 
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I. INTRODUCTION 


Clay minerals have a complex layered structure with exchangeable cations that can bind 
water molecules in the interlayers. With increasing pressure and temperature, or in the pres¬ 
ence of polar free radicals, these interlayer cations can be exchanged. This cation excahange 
capacity (CEC) of clay minerals affects their fluid conductivity,fIland permeabilitj^, dielectric 
permittivity.^ Complex dielectric properties of clay are crucial to determine hydrocarbon- 
contents in oil-rich rocks. These measurements have been made predominantly at frequencies 
in the kHz-rang^^ and between 0.5 MHz and 1.1 GHz.^ Complex conductivity of clayey 
materials between 1 milli-Hertz (mHz) and 45 kHz for CEC effects has been characterized 
and modeled.^ Canan^^ researched dielectric properties of smectite clay samples in detail, ex¬ 
plaining interlayer polarization and relaxation mechanisms between 30 kHz-300 MHz. Some 
clay minerals can swell due to hydration with water adsorbed in the interlayer depends on 
the charge of the interlayer cations. Electrical measurements can yield cation mobility. Di¬ 
electric measurements can be instrumental to characterize the water absorbed in smectite^ 
The conductivities of smectite clays saturated by monovalent cations has also been studied in 
detail.!® These frequency-and temperature-dependent measurements were in the frequency 
range between 300 Hz and 10 kHz, and from -150°C to -|-30°C. In a much higher frequency 
domain, the THz dielectric constants of layered silicates including muscovite, vermiculite, 
phlogopite, and biotite, have been measured by THz-time domain spectroscopy.!^ Unfortu¬ 
nately this technique is fairly noisy and low resolution since the THz pulses are very weak. 
So, in this paper we look at Monmorillionite using much higher resolution CW methods 
based on harmonic multiplication of phase stabilized microwaves, electronically generated; 
these methods provide very low noise/high dynamic range out to about 1.4 THz at present. 


II. METHODS 

At CSM we use three millimeter wave (or sub-THz) modalities to extract material prop¬ 
erties: (1) a quasi-optical System,!^!® to study bulk properties, (2) a near field scanning 
system,!^ to measure local properties, and (3) the open hemispherical cavity resonator,!^^^!! 
for samples that are too thin or too low-loss for quasi-optical techniques. In this work, using 
cavity resonance perturbation, we extract the complex dielectric constants of clay-thin Elms 
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in 100-165 GHz and investigate electrical properties in the presence of Ca’'"''/Na+-ions. We 
study how these cations influence conductivity of free carriers, and relaxations. We also 
compare our data with low frequency measurements.^ 

We use an open hemispherical cavity resonator with VNA (Vector Network Analyzer)^^ 
to measure electrical properties of thin sections of clay-samples with Ca+’''/Na'''-ions infused. 
The cavity is a structure with two copper mirrors positioned at certain distance (the “cavity 
length”) without any sidewalls. The top mirror is hemispherical and connected to two WR- 
10 waveguide couplers working as a transmitter and a receiver, and on the other hand the 
lower mirror is flat and smaller than the upper one in size. We measured the real part of 
refractive index of ~ 1 mm-thick glass substrate (borosillicate) to be 1.98 at 310 GHz which 
is the same as its published value.^ For details on the cavity and methodology see Rahman 
et al^. 

A. Open Cavity Resonator 

The principle of this technique is cavity perturbation. The changes in axissymetric mode 
prohles, mainly the frequency-shift and linewidth-variation, between an empty cavity mode 
and the same mode in presence of a sample, allows us to determine the complex dielectric 
constant of the sample.^^^^ The unloaded (empty) cavity has an axisymmetric mode spacing 
that is c/2L, where L is the distance between two mirrors, also known as cavity length. In 
our cavity, since L is around 15 cm, the unloaded mode spacing is about IGHz. 

Putting a sample on the bottom mirror perturbs the modes in a calculable but nontrivial 
way. To avoid geometrical factors, we do a second perturbation which involves flipping the 
sample upside down. Since the boundary values of the E-£eld are different, we are able to get 
a simple (geometry-free) formula for the complex permittivity from 3 sweeps around a OOq 
mode. This complex dielectric permittivity is related to total electrical/optical conductivity 
and absorption coefficienlP^ with the use of basic theory of electromagnetism as 

HreOi = 1207rcrre = SOueim 

where n^e and are are the real part of the refractive index and conductivity (in G'^cm”^), 
Sim is imaginary part of the complex dielectric constant, and uj and a are frequency and 
absorption coefficient expressed in terms wave numbers (cm“^ relative to the speed of light). 
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FIG. 1. Perturbations with substrate only, film up, and film down set-ups for clay with (a) Ca+"*“ 


and (b)Na^-ions (at room temperature) 


so lcm“^ corresponds to 30 GHz. The electrical conductivity is a macroscopic quantity which 
can be due to bound electrons (which is responsible for polarization and relaxation) or free 
electrons (which is Drude-type mechanism). This is also valid for absorption coefficient. 
These parameters, are and a, essentially describe the loss mechanisms in a material. 


B. Measurements 

By sweeping the VNA, we identified the axissymetric empty cavity modes based on 
constant frequency-spacing. This fixes the cavity-length, L, to be 145.56 mm and kept it 
unchanged throughout the experiments. Using a least-squares fit of the data to a Breit- 
Wigner model (as described in Rahman et we retrieve the eigenfrequency-shifts and 

modal quality factors (Q-values related to a linewidth) for substrate-only, film up and film 
down positions in order to apply the differential method. In order to determine uncertainty 
in the experiment, we repeat the entire procedure of inserting the sample, performing the 
measurements, and taking it out six times and calculating the variations in frequency shifts 
and linewidth changes of substrate. In Rahman et al^ we showed that by re-doing the 
entire procedure six times for borosillicate glass substrate, we obtained standard deviations 
less than 1.0% in the complex dielectric constant. This uncertainty also propagates to the 
calculations of dielectric constant. For consistency we confirmed during each trial that the 
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same part of the samples is probed. 

The frequency shifts due to him up, him down, and substrate only are used in the 
Eqs.([^ to determine the real part of the refractive index of a thin hlm.^^^*® This is called 
the diherential method since it involves measurements with the sample him up and down; 
the him will then see a diherent held because of the boundary conditions. The result is that 
sample geometry factors cancel out, greatly simplifying the calculation. Let nj and Ug refer, 
respectively to the refractive index of the him and substrate being measured. Then one can 
show that 

dUg — 1 

With 

SUf = l^{fup) ~ ^{s)- 
^(fdown) ^(s)' 

Here, iy{fup) , i^(fdown) , i^{s) represent the eigenfrequency associated with the him on the 
top (him up), him at the bottom (him down), and the substrate only, respectively. The 
term, 5vs, stands for the diherence between the eigenfrequencies associated with him at 
the bottom of the substrate and the substrate only, and the term, is the diherence 
between the eigenfrequencies with the him at the top of substrate and the substrate only. 
The required condition is the real part of the refractive index of substrate must be known. 
(Or since the substrate is large enough, we can measure its permittivity directly using 
the quasi-optical methodology. The real part of the dielectric constant of thin him can 
obtained by = We need Q-values of him up and substrate (only) to calculate the 

imaginary part of complex dielectric constant of the thin him. The Q-value of a resonant 
peak (perturbation) is dehned as Q=z/o/Ai/. The Q-value is related to the imaginary part of 
the refractive index by njm=l/2Q, so for a thin him, The latter expressions 

assume that Q is not too small (say less than 10). In a cavity such as ours, the Q could 
be 10^ — 10® depending on technical details of the waveguide coupling. In the diherential 
method, qisubstrate)^_ we are able to compute the imaginary part 

of the complex dielectric constant of the thin him by 



( 1 ) 

( 2 ) 

( 3 ) 
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FIG. 2. Experimental and simulated XRD patterns from oriented aggregate preparations from 
Ca'*'^ and Na^-ion stabilized clay samples used in EM analysis. 

C. Samples and Sample preparations 

We studied smectite clay minerals (SWy) from the Clay Minerals Society {\\www.clays.org\\) 
that were treated to yield homoionic, univalent (Na"*") and divalent (Ca^'*') samples. The 
samples were treated to remove carbonate and iron oxide cements with a Na-acetate buffer 
and Na-dithionite respectively.^ The sub-0.5 /xm equivalent spherical diameter size fractions 
were separated from bulk samples by standard centrifuge methods. Then the Na"*" saturated 
clay was thoroughly cleaned with dialysis to remove excess salt. To prepare the Ca^"'' 
exchanged smectite, a dialyzed Na"*" sample was treated with a IM solution of CaCh, and 
shaken for at least two hours, excess solution was decanted and the process was repeated 
twice more. The excess Ca^’*' salt was removed with dialysis.^ Oriented aggregates were 
made by evaporation onto glass slides to provide a sample ~ 4 cm long with at least 10 
mg/cm^ of clay.l221rhis evaporation method yields deposition of clay particles on the glass 
slide with the clay layers parallel to the glass slide.*^ 
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FIG. 3. Gravimetric water adsorption isotherms for SWy-Ca and SWy-Na samples following drying 
at 60G for 3 hours. 

XRD measurements showed systematic differences between Na"*" smectite and Ca’*"’'' smec¬ 
tite. The structure of smectite determined from oriented aggregate sample preparations and 
X-ray diffraction (XRD) prohle modelling and especially the interlayer water complex under 
different environmental conditions, is described in detail by Ferrage et q/JMIESIISSEZI^ Perrage 
et alW^ discuss the influence that the structure and organization of interlayer water has 
on the materials dielectric constant. For this reason, the diffraction prohles from the Na"*" 
and Ca’*"*' exchanged preparations used in the dielectric measurements in this study were 
simulated following the same methods as described by Drits and Sakharov*^ and Drits and 
Tchoubarl^. The experimental and modelled diffraction patterns for the Na"*" and Ca’''^-ion 
stabilized clay samples are shown in FIGj^ along with the thicknesses of bihydrated (2W), 
monohydrated (IW) and dehydrated (OW) layers and their proportions that were used in the 
simulated diffraction profiles. Under the ambient conditions of ~ 50% relative himidity (RH) 
in the laboratory the Na^ smectite is well described by random interstratihcations of 2W 
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(2%), IW (85%), and OW (13%). In contrast, the hydrated structure of Ca’''+ smectite was 
modelled by a different arrangement of layer thicknesses and proportions where 2W (81%), 
IW (9%), and OW (10%) were randomly interstratified. The XRD simulations in FIGj^ 
are not perfect with some misht present at different angular ranges for the Ca^’*' and Na"*" 
samples. This suggests that there is greater heterogeneity in the layer assemblage. However, 
the near perfect fit of the 001 reflection at 12.54 A (Na) and 15.05 A (Ca), and the generally 
good £t at higher angles indicate that the models represent the primary structures, with 
IW layers being dominant in the Na"*" sample and 2W layers being the dominant layer type 
in the Ca’*'^ sample. Smectite interlayer water as probed by XRD, considered crystalline 
water,^ may or may not all be coordinated to interlayer cations. Based on a combination 
of XRD and Grand-Canonical Monte Carlo simulationd^^ found that fluctuations of charge 
locations and water dipoles affect the dielectric constant. These authors also point out 
that gravimetric adsorption methods produce similar, but slightly higher molar equivalent 
water content than XRD experiments conducted under the same humidity conditions and 
based on the XRD simulation method. This is due to additional adsorption on surfaces and 
condensation in nm-scale pores. 


Thus, in addition to millimeter wave EM analysis, we performed a variety of measure¬ 
ments on the two samples, including thermal gravimetric analysis (TGA), gravimetric water 
adsorption, and subcritical nitrogen gas adsorption (SGA) to characterize physical differ¬ 
ences that exist between the Na"*" and Ca’*"’'' ex-infused samples. Na"*" smectite has higher 
SSA and dominant 1-3 nm pores typical for clay aggregates; in contrast, the Ga^"'' smectite 
sample has predominantly larger pores between 50-100 nm.^^SlAt ~ 50 RH gravimetric water 
adsorption data (FIG. from SWy-Ga^"'' and SWy-Na"*" showed that the Ga^"'' smectite 
adsorbed ~ 17.5 wt.% following dehydration at 60°G for 3 hours, compared to 10 wt.% for 
the Na"*" smectite under the same conditions. The hydration enthalpy is about four times 
greater for Ga^"*" than Na"*" cations and is consistent with these gravimetric data. However, 
in the oriented aggregate preparations used for the EM analysis the contrasting structural 
arrangement of the interlayer crystalline water between the Na"*" and Ga"*"’*' smectite samples 
(FIG. must play a significant role in the EM response.!^ 



FIG. 4. High frequency (100-165 GHz) data of (upper left) real part and (upper right) imaginary 
part of the complex dielectric constants of clay samples with Na"’“/Ca^'’“-ions (at room tempera¬ 
ture); (lower middle) conductivity of both samples (inner plot is the expansion of the lOO-llOGHz 
responses). 


III. RESULTS AND DISCUSSION 


Referring to FIGj^ one can see that the real part of the dielectric constants of both 
clay samples are nearly dispersion free (i.e., frequency independent) in our range of mea¬ 
surements, indicating that carrier concentrations are low in both.^^^® The imaginary part 
of dielectric constant and electrical conductivity of Na+-ionized sample, increase at two dif¬ 
ferent rates: faster in the higher frequency range of our measurement and slower at low 
frequencies. The presence of both Debye relaxation!^, and Maxwell-Wagner relaxatioiP^ 
at radio frequencies has been reported.® But between 100 and 165 GHz, phonon induced 
relaxation is dominant.^ 

Revil et a/.® proposed a theoretical model to study complex conductivity-dependences on 
cation exchange capacity (GEG), specific surface area (SSA), and salinity for clay samples at 
low frequency. This model also relates GEG to SSA with consistency.® It is reported that SSA 
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increases exponentially for Na’''-clay and linearly for Ca’''’''-clay sample.^ From our research, 
it is evident that conductivity depends on CEC or SSA. The larger CEC or expnonetial- 
growing SSA can contribute to more disperse Stern layer and the smaller CEC or linearly- 
progressive SSA stabilize the Stern layer. Thus, electrical conductivity is linked to CEC and 
zeta potential in the Stern layer.® The Na"*" makes a thicker unstable double layer where these 
high mobility ions are able to polarize rapidly. Therefore, the conductivity is more dispersive. 
This interlayer polarization is correlated to relaxation mechanisms. The relaxation processes 
are, therefore, dependent of ^-potential which is also correlated to CEC.^On the other hand. 
The imaginary part of dielectric constant and conductivity of the sample with Ca’''^ increase 
monotonically and sublinearly. The Ca^"*" creates more stable double Stern layer. Due to low 
mobility, the interlayer polarization is less disperse so its conductivity is sublinear. Since, 

= 2 {1/v) Urei the effect of (1/zz) is more into imaginary part of complex dielectric constant 
of samples with Na’*"’*' than that of Ca"*"^. In this case, at high frequencies faced steeper 
decrease than at lower frequency ends. 

Comparing the results from Raythatha and Sen® with ours (FIG§, it is clear that in 
the sub-THz, the electrical conductivities of clay-Na’''/Ca’'"''-ions are almost 2 to 3 orders 
of magnitude greater than in the RF/microwave range. On the other hand, the real part 
of dielectric constants for both samples, are roughly an order of magnitude larger in the 
RF/microwave range than in the sub-THz. The huge RF/microwave permittivities have been 
attributed to polarization of the double-layer surrounding the ions.® The lower (and nearly 
frequency-independent) values of e^e suggest that there are depletions of mobile charges due 
to phonon interactions causing the diffused layers to thin. 


IV. CONCLUSION 

We measured sub-THz complex dielectric properties of clay samples with Na’''/Ca'''’'', and 
compute their e^e and cr^e- We illustrate the connections between electromagnetic parameters 
and components of surface chemistry such as CEC or SSA, (^-potential. This also enables 
us to study the clay content and free radicals in shales, and to investigate their CEC and 
(■-potential-dependences. In future, we will study more different ionized clay samples and 
broader frequency ranges to capture more complete frequency-dependences of e^e and eim- 
This will allow us to model computationally the dispersions of these parameters. 
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